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Discovery of a massive variable star with Z=Z /36 in the 
galaxy DDO 68 



oo 
O 

o 

(N 

>y 
c3 



(N 
(N 

Si 

6 



in 
in 
cn 

in 
o 
oo 
O 



S. A. Pustilnik, 1 * A. L. Tepliakova, 1 * A. Y. Kniazev, 2 ' 1 * A. N. Burenkov 1 * 

1 Special Astrophysical Observatory of HAS, Nizhnij Arkhyz, Karachai-Circassia 369167, Russia 

2 South African Astronomical Observatory, Observatory Road, Cape Town, South Africa 



X 



Accepted 2008 May 1. Received 2008 April 26 



ABSTRACT 

The Local Volume dwarf galaxy DDO 68, from the spectroscopy of its two 
brightest HII regions (Knots 1 and 2) was designated as the second most metal- 
poor star- forming galaxy [12+log(0/H)=7.14]. In the repeated spectral observations 
in 2008 January with the 6-m telescope (BTA) of the HII region Knot 3 [having 
12+log(O/H)=7.10±0.06], we find a strong evidence of a transient event related to a 
massive star evolution. From the follow-up observation with the higher spectral res- 
olution in 2008 February, we confirm this phenomenon, and give parameters of its 
emission-line spectrum comprising of Balmer HI and He I lines. The luminosities of 
the strongest transient lines (Ha, H/3) are of a few 10 36 erg s _1 . We also detected 
an additional continuum component in the new spectrum of Knot 3, which displays 
the spectral energy distribution raising to ultraviolet. The estimate of the flux of this 
continuum leads us to its absolute V-band magnitude of ~-7.1. Based on the spectral 
properties of this transient component, we suggest that it is related to an evolved 
massive star of luminous blue variable type with Z=Z Q /36. We briefly discuss ob- 
servational constraints on parameters of this unique (in the aspect of the record low 
metallicity of the progenitor massive star) event and propose several lines of its study. 
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individual:DDO 68 (UGC 5340) - stars: emission-line, Be - stars: variables: other 



1 INTRODUCTION 

The most metal-poor galaxies [or extremely metal- 
deficient, XMD, wjth 12+log(0/H)<7.65, e.g., review by 
iKunth fc OstlirJ (|200Ch ] are considered as the best lab- 
oratories to study in great detail the processes related 
to galaxy evolution and star formation (SF) in young 
galaxies at high redshifts. Among several XMD galax- 
ies with O/H near the bottom of dwarf galaxy distri- 
bution, namely those with parameter 12+log(0/H)<7.3, 
DD O 68 (UGC 5340) with its value of 7. 14 is the nearest 
one (IPustilnik. Kniazev. fc Pramskiil[2005l : llzotov fc Thuanl 
120071 . hereafter PKP and IT07). Its distance was accepted 
as 6.5 Mpc in PKP. However, for its Vh c n o =502 kms" 1 , 
due to the large ne gative pecu l iar v elocity correction 
in this space region l|Tullv et al.l 120081 ). this can be as 
large as ~10 Mpc. Nevertheless, this galaxy is one of 
the most perspective to study individual massive stars 
with the future Extremely Large Telescopes (ELTs) (e.g., 
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iKniazev fc Pustilnikl [20061 ) . In the recent papers DDO 68 
was sho wn to be a likely result of a merger of two gas-rich 
dwarfs l|Ekta. Chengalur. fc Pustilnikl r200sT) and to show 
up colours on the SDSS images, which for the standard 
Salpeter IMF are consistent with no tracers of stars older 
than ~1 Gyr (jPustilnik. Tepliakova, fc Kniazevll2008l ). The 
current O/H value of DDO 68 is based on the direct mea- 
surements by the classic T c -method of two brightest HII 
regions in the Northern ring, called Knots 1 and 2 (PKP, 
IT07). In this Letter we present new BTA spectra of Knots 
1 and 3. The latter displays transient emission lines and blue 
continuum, which we attribute to a massive variable star of 
luminous blue variable (LBV) type. 



2 OBSERVATIONS AND DATA REDUCTION 

The long-slit spectral observations of Knots 1 and 3 (15" 
in between) were conducted with the multimode instru- 
ment SCORPIO l|Afanasiev fc Moiseevl 120051 ) installed in 
the prime focus of the SAO 6 m telescope (BTA) on the 
nights of 2008 January 11 and February 4. The grism 
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VPHG550G was used with the 2Kx2K CCD detector 
EEV 42-40 on 2008 January 11 and the grism VPHG1200G 
on 2008 February 4. These set-ups gave the range 3500- 
7500 A with ~2.1 A pixel -1 and the full width at half- 
maximum (FWHM) of ~12 A and the range 3900-5700 A 
with ~0.9 A pixel -1 and FWHM of ~6 A along the disper- 
sion, respectively. The scale along the slit (after binning) was 
0'.'36 pixel -1 . Six and five 15-min. exposures were obtained 
on nights 2008 January 11 and 2008 February 4, respec- 
tively, under the seeing of 1.2 and 1.3 arcsec. The objects 
spectra were complemented before or after by the reference 
spectra of He-Ne-Ar lamp for the wavelength calibration. 
Bias and flat-field images were also acquired to perform the 
standard reduction of 2D spectra. Spectral standard star 
Feige 34 (jBohlinl Il996h was observed during the night for 
the flux calibration. 

All data reduction for Knots 1 and 3 and emission line 
measurements for Knot 1 were performed similar to that de- 
scribed in PKP. Namely the standard pipeline with the use 
oflRAFQand MIDA93 was applied for the reduction of long- 
slit spectra, which included the next steps: cosmic ray hits 
removal, bias subtraction, flat-field correction, wavelength 
calibration, night-sky background subtraction. Then, using 
the data on spectrophotometry standard stars, all spectra 
were transformed to absolute fluxes. ID spectra were ex- 
tracted by summing up, without weighting, of 12 and 6 rows 
along the slit, in Knots 1 and 3, respectively. The emission 
lines with t heir errors were meas ured in the way described 
in details in iKniazev et ail l|2004h . 



3 RESULTS 

In Fig. [T] the relevant new and old spectra, obtained with 
the same set-up, are shown. The bottom panel displays 
the spectrum of Knot 3 obtained on 2008 January 11 with 
FWHM'~12 A, while the middle panel shows the similar 
spectrum of Knot 3 obtained on 2005 January 13. Signifi- 
cant changes in relative line intensities have occurred during 
the past ~3 years. In the top panel, we show two spectra of 
Knot 1, obtained on 2008 January 11 (upper) and 2005 Jan- 
uary 8. The latter spectrum is shifted down by 0.05 flux unit 
to allow better visual comparison of both spectra. They have 
very similar parameters of both relative line intensities and 
underlying continuum and show no evidences for systematic 
difference. 

In Table [TJ we present the line intensities of all relevant 
emission lines measured in the spectrum of Knot 1 of 2008 
January 11, normalized by the intensity of H/9 and corrected 
for the foreground extinction C(H/3) and equivalent widths 
of underlying Balmer absorption lines (EW(abs)), given in 
the bottom part of this table. 

We also present in this table the electron temperatures 
for zones of emission of [O ill], [O n] and [S ill], N e , and ionic 



1 IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatory, which is oper- 
ated by the Association of Universities for Research in Astron- 
omy, In. (AURA) under cooperative agreement with the National 
Science Foundation (NSF). 

2 MIDAS is an acronym for the European Southern Observatory 
package - Munich Image Data Analysis System. 



Table 1. Line intensities and derived parameters of Knot 1 
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C(H/3) dex 
EW(abs) A 
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S+/H+(xl0 -7 ) 
S++/H+(xl0 -7 ) 
ICF(S) 
log(S/0) 



0.00±0.03 
2.50±0.14 
39.1±0.7 
162±3.9 
476±33 

18.81±1.12 
15.31±0.86 
17.32±0.93 
378±294 

0.301±0.044 
1.192±0.167 
0.039±0.006 
1.532±0.173 
7.19±0.05 

0.201±0.027 

1.285 
-0.77±0.08 

0.527±0.058 
3.132±0.742 

1.601 
-1.56±0.13 



a in units of 10 16 ergs s 'cm 2 . 



abundances of oxygen, neon and sulphur (with respective 
ICFs - ionization correction factors), along with the total 
abundances derived for the above measured line intensities 
with the cla ssic T c method as described in IKniazev et alj 
(2004. 2005). The latter are consistent with the results from 
both PKP and IT07. The relative line intensities of Knot 1, 
measured on the spectrum 2008 February 4 (not shown here 
due to the lack of space) with a higher resolution but with a 
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more narrow range, are consistent with those from Table [T] 
However, they can not be used for the independent deriva- 
tion of O/H since this spectrum does not include the line [O 
n]A3727. 

One of the aims of new observations of Knot 3 was an 
attempt to improve its element abundances determination. 
However, due to appearance in the spectrum of 20008 Jan- 
uary 11 of additional 'broad' components and the difficulty 
due to low spectral resolution (FWHM~12 A) to disentan- 
gle nebular emission lines from variable 'broad' components, 
we should postpone this task till getting the spectrum with 
a higher resolution. The same relates to the spectrum of 
2008 February 4, with FWHM~6 A. This appears also in- 
sufficient to get reliable line intensities of nebular emission 
lines. Before to return to analysis of lines of the variable 
component, we briefly describe one more estimate of O/H 
for Knot 3, derived from its old spectrum, which will be used 
in Discussion. 

From the line intensities of Knot 3 in table 3 and fig. 2 
of PKP (also shown in the middle panel of Fig. [1} , with the 
use of direct T e method for the observation with the better 
S-to-N ratio in I([0 IIIA4363]) (2005 January 13), the derived 
value of 12+log(0/H), assumed to be related to the metal- 
licity of the variable star, discussed below, was 7.11±0.11. 
This motivates an attempt to improve the accuracy of O/H 
in Knot 3. 

We apply below the so-called 'semi-empirical' method, 
suggested by IT07. It is intended for abundance determi- 
nation in the case when the principal faint line [Oiii]A4363 
is too noisy or undetected. This uses the rather tight cor- 
relation between T(0 in) and the total intensity of the 
strong oxygen lines relative to I(H/3) (Stasihska & Izotov 
2003). After T(Oiii) is estimated from the total intensity 
of [On]A3727, [O ill] AA4959, 5007 lines via the relation pre- 
sented in IT07, all other calculations are made as in the clas- 
sic T e method, using measured line intensities. The robust- 
ness of the semi-empirical method in the very low-metallicity 
regime was tested in IT07 and by us, using all available data 
on the DDO 68 HII regions. These O/H estimates relative to 
those obtained by the direct T e method, appeared very sta- 
ble and showed a scatter consistent with the errors of the di- 
rect method determination and those for the semi-empirical 
method itself. Applying this method to the same data on line 
intensities, and ignoring the faint line [O ill] A4363, we obtain 
for Knot 3 the value of 12+log(O/H)=7.10±0.06, which is 
nicely consistent with the value derived from the direct T e 
method. 

We therefore accept for further discussi on this value of 
O/H , which corresponds to 1/36 of (O/H)© l|Asplund et alj 
12004 ) . As the nearest of young star-forming regions (age of 
5.0 Myr, PKP) with such a low metallicity, Knot 3 is cer- 
tainly one of the most attractive targets for the future Ex- 
tremely Large Telescopes (ELTs) to study individual mas- 
sive stars with metallicities typical of very young galaxies in 
the early Universe. 

Fortunately, by good luck, some tracers of the indi- 
vidual massive stars in nearby metal-poor galaxies can be 
found occasionally in the spect ra of star-forming regions 
(|lzotov. Thuan. fc Guseval 120071 ). This allows to conduct 
their studies with modern large telescopes. DDO 68 appears 
to be one of such galaxies, but outstanding for its the record 
lowest metallicity. In Fig.[TJ we display the old and new spec- 



Table 2. Parameters of residual lines in Knot 3 
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tra of Knot 3, with resolution of FWHM~12 A, in which 
pronounced differences have occurred during the past 3 yr. 
To better describe the variable component in the spectrum 
of Knot 3, including its continuum, we show its difference 
spectrum in the top panel of Fig. [2] The normalization of 
the old spectrum (by a factor of 0.85) was performed before 
the subtraction, in order to fully remove the (assumed) HII 
region nebular lines [O iii]A4959,5007. 

In the bottom panel of Fig. [2] we show the spectrum 
of Knot 3 with resolution of FWHM~6 A, in which broad 
components of Balmer and He I lines are more evident. The 
subtraction of nebular lines (shown in red) in this spectrum 
was performed by scaling the same lines in the spectrum 
of Knot 1, exposed on the same long slit, just ~15 arcsec 
a side. We present the parameters of 'narrow' and 'broad' 
components of the residual lines, shown in black, in Table 
where the names of the columns are self-explaining. To 
determine parameters of 'narrow' and 'broad' components 
(index n and b for Balmer and Hel lines) we performed the 
procedure of two-component Gaussian fitting in MIDAS. To 
understand the relative velocities of the transient emission 
and that of Knot 3 HII region, the heliocentric velocity of 
the latter, defined on the spectrum with FWHM=6 A, is 
V(Knot 3)=564±27 kms" 1 . 



4 DISCUSSION 

The natural interpretation of this transient emission-line 
component is the appearance of a luminous variable star. 
Some of the observed and derived parameters of this unique 
variable star and the related star-forming region are as fol- 
lows. The integrated l/-band magnitude of Knot 3, as mea- 
sured on the respective image in PKP, is ~21.2. The latter 
is consistent with V calculated through convolving the 2008 
January spectrum of this region with the ^/-filter bandpass, 
namely V=20.94. A similar convolution of the residual spec- 
trum in the top panel of Fig. [2] allows us to estimate the 
^/-magnitude of the transient to be 22.95. Its continuum, 
clearly rising to the ultraviolet (UV), indicates a hot star 
of OB spectral type. For the accepted distance module of 
DDO 68 /x=30.0 (D=10 Mpc) (see Pustilnik et al. 2008), 
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Figure 1. Top panel: The spectra of DDO 68 Knot 1 on 2005 
January 8 and 2008 January 11, superimposed, with the former 
one shifted down in flux by 0.05 unit. The spectra look very sim- 
ilar, including both lines and continuum. Middle panel: Spec- 
trum of Knot 3 in DDO 68 on 2005 January 13. Bottom panel: 
Spectrum of Knot 3 in DDO 68 on 2008 January 11. Pay at- 
tention on the significant enhancement in the latter spectrum of 
He I lines (A4471, A4713, A5876, A6678 and A7065), the appear- 
ance of asymmetric broad wings in the strongest Hel lines, and 
the change of relative intensities of [O III] A4959, 5007 and Balmer 
lines. 



the absolute magnitudes, Mv, of Knot 3 and the variable 
star are -9.1 and -7.1, respectively. 

For the estimated Knot 3 age of ~5 Myr (PKP), one 
can calculate the mass of its young stellar cluster, using, 
e.g., models for evolving stellar populations of PEGASE2 
package (Fioc & Rocca-Volmerange 1999). Accounting for 
~20 per cent light contribution of nebular emission and a 
small correction for the Galaxy extinction, with the accepted 
stellar metallicity of z=0.0004 and the standard Salpeter 



Figure 2. Top panel: The difference spectrum of Knot 3 with 
resolution of ~12 A, obtained by subtraction of the two spectra 
taken 3 years apart. The residual features represent the massive 
star wind emission lines, while the underlying continuum presum- 
ably indicates the spectral energy distribution of a variable hot 
star. The marked lines are H8+He I 3889, He 3970, He I 4026, H<5 
4102, H7 4340, He I 4388, He I 4471, He I 4713, H/3 4861, He I 4922, 
He I 5016 (affected by subtraction of [Oin]5007), Hel 5876, H«, 
Hel 6678. Bottom panel: The spectrum of Knot 3 with resolu- 
tion of ~6 A. Shown here are both the original spectrum (red) 
and a difference spectrum (black), obtained by the subtraction 
of scaled lines in the spectrum of Knot 1, registered on the same 
slit. Continuum remains original. Broad wings of both Balmer and 
He I lines are visible as well as an indication of P Cygni profiles 
in Balmer lines. 



Initial Mass Function with Mi ow and M up of 0.1 and 120 
M , we obtain M c i ustcr = 8.5xl0 3 M . 

To further address the type of a massive star related to 
the observed transient, we briefly describe its emission-line 
properties. In both spectra of the variable massive star in 
Fig. [2] only the Hydrogen Balmer series and He I emission 
lines are detected. This is in agreement with its extremely 
low metallicity of Zq/36. The observed line profiles, with 
broad asymmetric wings and an indication of P Cygni profile 
resemble those of B-supergiant LBVs. These represent the 
short transition phase in the post- main- sequence evolution 
of massive stars with M<50 Mq. The latter are exemplified 
by spectra of stars He3-519, AG Car, 77 Car and others in 
Kudritzki (1998) and Walborn & Fitzpatrick (2000). Their 
'narrow' components are close in radial velocity to that 
of the related HII region. However, their widths (FWHM) 
of ~700 kms^ 1 are clearly larger than the instrumental 
ones. The nature of both 'narrow' and broad components 
(FWHM~1000-2000 kms" 1 ) is not clear. For the accepted 
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DDO 68 distance of 10 Mpc, the H/3 line fluxes correspond 
to luminosities of ~4.6 and ~0.8 10 36 erg s _1 in 'narrow' and 
broad componen ts, respectively. Both elec tron scattering in 
dense envelopes (jBernat fc Lambertll 19781 ) and stellar winds 
could contribute to their appearance. The asymmetry of the 
broad components can be partly due to absorption of the 
approaching region of the (non-spherical) wind. This has im- 
portant implications, since in very low-metallicity stars the 
line-driven wind can be inefficient, and continuum-drive n 
eruptions might be more significant (|Smith fc Owocki 2006). 
The latter, in turn, will have important implications for evo- 
lution of Population III stars and the early metal enrichment 
of the Universe. 

The hunting for the most metal-poor local galaxies as 
useful testbeds for models of young galaxies was rather suc- 
cessful during the last years, with the lowest HH-region oxy- 
gen abundance found so far in SBS 0335-052 W on the 
level of 12+log(O/H)=7.12±0.03 (llzotov, Thuan. fc Guseval 
2005). The star-forming region Knot 3 of a much closer 
galaxy DDO 68 appears to have very similar metallicity of 
12+log(O/H)=7.10±0.06. DDO 68 displays also other un- 
usual properties, such as HI and optical morphology and 
kinematics of a likely recent merger and a rather young 
{SX Gyr) the oldest detectable stellar population (PKP), 
|Ekta et al.ll 20081 ; IPustilnik et al.ll2008r ). 

Due to its relative proximity and very low metallicity, 
this galaxy and its young SF regions are among the best 
targets for detailed studies of star- formation, its feedback, 
and massive star properties typical of high-redshift young 
objects. Monitoring of its several young star-forming re- 
gions, with ages of 4 to 7 Myr, to discover variable mas- 
sive stars similar to the one reported here, will be very 
important. The subsequent spectral studies of such ob- 
jects, could lead to interesting findings in extremely metal- 
poor environments and the properties of massi ve metal-poor 
stars , actively studied in recent year s (e.g., iMevnet et al.l 
l2007i ; lYoon." Cantiello. fc Langerll2008l ). In particular, imag- 
ing with the HST of Knot 3 would help to localise this lumi- 
nous star (y~23), and provide its spectral energy distribu- 
tion in broad range from UV to NIR. This information will 
be helpful for spectral study of this star with the next gen- 
eration space telescope JWST. Even its medium resolution 
spectra (R~5000) would provide the important constraints 
on rotation of such stars, while NIR spectra might shed light 
on the properties of related wind/shells. 

Massive low-metallicity stars are the main suppliers of 
energy, momentum and metals duri ng the early evolution o f 
galaxies in the young Universe fe.g-. lBarkana fc Loebll200lf ). 
Their properties and details of their evolution are important 
for estimating the effect of SF feedback on galaxy evolution 
and their effect on re-ionization of the Universe. Thus, the 
discovery and the study of such metal-poor massive stars 
with the next generation of ELTs is directly relevant to in- 
sights in the early Universe evolution. 

Summarising the results and discussion above we draw 
the following conclusions: 

(i) The spectrum of HH-region Knot 3 in the nearby XMD 
galaxy DDO 68 changed substantially during the period be- 
tween 2005 January and 2008 January, displaying enhanced 
intensities and broad components of Balmer and He I lines, 



as well as an additional continuum component typical of a 
hot star. 

(ii) The luminosity, blue continuum and emission-line pa- 
rameters of this variable emission resemble those expected 
for a LBV with metallicity of Z=Z0/36 (as measured on 
nebular emission in Knot 3). This suggests the discovery of 
the first individual massive star with the record low metal- 
licity and implies good prospects for study such stars with 
the future ELTs. 
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